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Three-Loop Mixing of Dipole Operators 
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We calculate the complete three-loop 0{a 3 a ) anomalous dimension matrix for the dimension- five 
dipole operators that arise in the Standard Model after integrating out the top quark and the heavy 
electroweak bosons. Our computation completes the three-loop anomalous dimension matrix of 
operators that govern low-energy |A_F| = 1 flavor-changing processes, and represents an important 
ingredient of the next-to-next-to- leading order QCD analysis of the B — > X a ^ decay. 
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Weak interaction phenomena at energies much below 
the electroweak scale are most conveniently described 
in the framework of an effective theory that is derived 
from the Standard Model (SM) by integrating out the 
top quark and the heavy electroweak bosons. The La- 
grangian of such an effective theory 



£eff = -CqCDxQED + £weak , 



(1) 



is a sum of the conventional QCDxQED Lagrangian for 
the remaining SM fields and a linear combination 



^-■wcak 



(2) 



of dimension > 5 operators Qi that are built out of those 
light fields. Ci(fi) are the corresponding Wilson coeffi- 
cients that depend on the renormalization scale fx. 

For most phenomenological applications, only opera- 
tors of dimensions five and six are relevant. The com- 
plete set of such operators consists of: %) dimension-six 
four-fermion operators, ii) dimension-six purely gluonic 
operators, and Hi) dimension-five dipole operators 



and 



V^TVG" . (3) 



Here u^ v = i/2 [7^, 7"], while tp and ip' stand for fermion 
fields of opposite chiralities. Their flavors may or may 
not be the same. The electromagnetic and strong field 
strength tensors are denoted by and G^, respec- 
tively. T a are the SU(3)c generators for the consid- 
ered fcrmions. For off-shell calculations, additional oper- 
ators that vanish by the QCD x QED Equations of Motion 
(EOM) must be included (see below). 

The structure of £ wca k remains the same in any 
SU(3)c x U(l) em gauge-invariant extension of the SM 
that does not contain exotic light bosons. New physics 
effects are therefore encoded in the values of the Wilson 
coefficients only. 

The dipole operators introduced in Eq. © are relevant 
in a variety of phenomenological applications, ranging 
from electric and magnetic moments of the leptons and 
nucleons to radiative decays, such as [i — > ej, ft — > E7, 
B -> K*y, B -> frf and B -> X s j. 



The Wilson coefficients are determined by matching 
Green's functions of the effective theory and the SM (or 
its extension) at the electroweak (or higher) scale /ihigh- 
Next, one applies the Renormalization Group Equations 
(RGE) 



~jd/')<"./(/') 



(4) 



to evolve Ci(/i) to the relevant low-energy scale /xi ow - In 
this way, large logarithms ln(A*high/Miow) are resummed 
from all orders of the perturbation series. 

Neglecting QED effects, the Anomalous Dimension 
Matrix (ADM) j(fi) has the following perturbative ex- 
pansion 
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fc=0 



4.7T 



k+1 



7 



CO 



(5) 



where a s = <7 2 /(47r) is the strong coupling constant. 

The purpose of our present work is to evaluate the en- 
tries of 7^ that correspond to the 0(a 3 ) strong mixing 
of the dipole operators of Eq. containing quark fields. 
These entries can be extracted from the three-loop QCD 
renormalization constants in the effective theory. Lower- 
order entries are already known from previous calcula- 
tions 0, 0, 0, II 0. The three-loop QCD self-mixing 
of the electromagnetic dipole operator coincides with the 
anomalous dimension of the rank-two antisymmetric ten- 
sor current that has been calculated in Ref. Q . We con- 
firm all these findings. Our remaining three- loop results 
are entirely new. 

The main phenomenological motivation for our work 
is to provide a new contribution to the calculation of 
the B — ► X s j branching ratio at the Next-to-Next-to- 
Leading Order (NNLO) in QCD. Including these 0(a 2 s ) 
corrections is necessary to reduce the theoretical uncer- 
tainty of the SM calculation Q below the current ex- 
perimental one ||- Several steps in this direction have 
already been made 0, 0, [n| . Our findings can also be 
relevant for the CP-odd electric dipole moment of the 
neutron, provided the new physics matching scale is suf- 
ficiently high. 
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Our calculation completes the three-loop QCD ADM 
for the whole dimension-five part of the effective La- 
grangian, because the dipole operators in Eq. © are the 
only EOM-non- vanishing operators in this sector. Simul- 
taneously, our results establish the whole three-loop QCD 
ADM for the \AF\ = 1 operators of dimensions five and 
six that arise in the SM case — all other three-loop ADM 
entries for such operators are known from previous pub- 
lications 0, 0, 0| . The only ADM entries that remain 
to be calculated for B -> X s -/ at the NNLO in QCD 
correspond to the four-loop mixing of certain four-quark 
operators into the dipole operators . 

For dcfinitcness in the further discussion, we shall 
choose the flavors, chiralitics, normalization and names 
of the dipole operators as it is usually done in the phe- 
nomenological analyses of B — > X s ^ 1 namely 



2 -m b (s L a^T a b R ) G% . 
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However, we stress that the 2x2 ADM which we calculate 
is the same for any pair of such quark dipole operators, 
including the flavor-conserving ones. There is no mixing 
between dipole operators of different flavor content, even 
in the flavor-conserving sector 0, |(| . 

In order to remove the divergences of all possible 
off-shell one-particlc-irrcducible (1PI) Green's functions 
with single insertions of Qj and Qs, we have toj^troducc 
the following EOM-vanishing countcrterms 



Qjpip 



1 



16tt 2 



mbS L ppb R . 



Qp0 = Y^ m b^ (j><p-<pp)b R , 



(7) 



where D u = d u + ig s G u + ieQ d A u and D u = d u -igsG^ 
— ieQdAfj, denote the covariant derivatives of the gauge 
group SU(3)c x J7(l) cm acting on the fields to the right 
and left, respectively, and we have used the definition 
= G^T a for the matrix-valued gluon field. A u is the 
photon field, and the color generators are normalized so 
that Tr(T a T h ) = S ab /2. 

Notice that the operator Qjp]pis gauge- invariant, while 
Qip0 is not. The appearance of such operators is ex- 
pected on general grounds [l4l fl5| . In principle, one 
could also encounter nonphysical countcrterms that can 
be written as Becchi-Rouct-Stora-Tyutin (BRST) varia- 
tions of some other operators, so-called BRST-exact op- 
erators. However, they turn out to be unnecessary in the 
case of the dipole operator mixing. This issue is discussed 
in more detail in Refs. |3.l5LHl|. 

We perform the calculation using dimensional regular- 
ization and the MS scheme. As far as the matrix 75 
is concerned, its only relevant property in our case is 
[757 7^7^] = 0, which, to our knowledge, holds in all the 
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FIG. 1: Typical examples of 1PI diagrams describing the self- 
mixing of Q7 at the one-, two- and three- loop level. 



commonly used schemes for the treatment of 75, includ- 
ing the Naive Dimensional Regularization (NDR) and 
t'Hooft-Vcltman (HV) schemes. However, the ADM be- 
yond one loop in the Dimensional REDuction (DRED) 
scheme is different from the one we find here because this 
scheme does not coincide with standard dimensional reg- 
ularization even in the absence of 75. A description of 
the properties of the NDR, HV and DRED schemes, as 
well as a list of relevant original articles can be found in 
Ref. [H. 

The necessary three-loop renormalization matrix is 
found by calculating the one- and two-loop b — > s, 
6 — » 57, b — > sg amputated Green's functions with sin- 
gle insertions of Q7, Qs, Qipjp and Qipg, as well as the 
three-loop b — * 57 and b — * sg amplitudes with insertions 
of Qr and Qs. Sample diagrams are shown in FIG. ^ 
The corresponding one-, two- and three- loop amplitudes 
are evaluated using the method that has been described 
in Refs. 0, 0- E3- We perform the calculation off shell 
in an arbitrary gauge, which allows us to explicitly 
check the gauge-parameter independence of the mixing 
among physical operators. To distinguish between in- 
frared and ultraviolet (UV) divergences, we introduce a 
common mass M for all fields, expanding all loop inte- 
grals in inverse powers of M. This makes the calcula- 
tion of the UV divergences possible at three loops, as M 
becomes the only relevant internal scale, and three-loop 
tadpole integrals with a single nonzero mass are known 
0,^3- On the other hand, this procedure requires to 
take into account insertions of the nonphysical operators 
Qjpip and Q jpg , as well as of the following counterterm 
of dimension three: 



M 2 m b s L b r 



(8) 



A comprehensive discussion of the technical details of 
the renormalization of the effective theory and the actual 
calculation of the operator mixing is given in Refs. 0,0. 

Having summarized our method, we now present our 
results for arbitrary numbers of down- and up-type quark 
flavors denoted by rid and n u , respectively. The ADM 
depends on the total number of active quark flavors nj = 
nu + nd, and their "total" electric charge Q = n u Q u + 
n-dQd- The regularization- and renormalization-schcmc 
independent matrix 7^°) is given by 
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While the matrix •yW is renormalization-scheme inde- 
pendent, j' 1 ' and 7^ are not. In the MS scheme, we 
obtain 



7 (i) 



1936 224 „ 
9 27 / 







( 368 224 „ \n . 1456 61 r 
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(10) 



and 

(2) _ 



7 



Q_(l856 + 1280„ / ) Qd 

352 „2 
"ST"/ 

352 2 \n , 



1 _ 1312 



C(3) 



17108 _ 



(11) 



We remark that the explicit electric charge Qd originates 
solely from the quarks in the operators, and thus has 
to be replaced by Q u for operators containing up-quark 
fields. As it is characteristic for three-loop anomalous 
dimensions, the entries of 7^ involve terms proportional 
to the Riemann zeta function £(3) ps 1.20206. 

Of course, the presence of the bottom quark mass in 
the normalization of the dipole operators of Eq. JBJ af- 
fects the values of 7'*). Had we decided to define the 
operators without quark mass in their normalization, the 
results in Eqs. Q to (fTTJl would need to be replaced by 

( fc ) ( k )x u 
7h - 7™ dij, where 
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280 



— C(3) ) n, - — 



are the expansion coefficients of the quark mass anoma- 
lous dimension. In particular, verifying that our result 
for yffi is in agreement with Eq. (8) of Ref. requires 
to perform such a replacement because no quark mass 
was present in the normalization of the tensor current 
considered there. We note that there is a misprint in the 
last line of Eq. (7) of the latter paper. Obviously, the 
factor Cp should read Nj. 

The RGE for the Wilson coefficients given in Eq. |@J, 
has the following general solution 



Ci(nb) = Uijbib, fi w )Cj(fi w ) 



(13) 



where the matching and the low-energy scales have been 
denoted by fi w and fib, respectively. In the B — » X s y 
case, one has fi w = 0{M W ) and m, = 0(m&). The evo- 
lution matrix Uij(fib, fiw) depends on the strong gauge 
coupling ratio 77 = a s (fi w )/a s (fib). The ADM that we 
have calculated allows us to find the complete 0{a 2 s ) con- 
tributions to 

AC 7 {nb) = fyjiVb^wWjifiw), (14) 
1=7,8 



and 



ACs(nb) = Ussinb, fi w )C 8 {fiw) 



(15) 



Denoting these 0{ot 2 s ) contributions by A' 2 ' C7 (fib) and 
A( 2 'Cs(^&), respectively, and u sing the general NNLO 
formalism presented in Rcfs. |llHl9j . we obtain for n/ = 
5, Q d = -1/3 and Q = 1/3: 



37208 / 39 62 \ (i). , 
~ "476T V ~ 7,23 J 7 0v) 



#4 2) (mw) 



/ 7164416 
V 357075 

256868 a. 

-77 2: 



297664 39 
7723 



14283 
6698884 



(16) 



n* KVw) 



14283 ' 357075 

I6.65I677M - 61. 0768 772! + 44.4252 77^) Cf\fi w ) 
+ (^36.4636 77^ -44.4043 7^ - 135.3141 77^3 

+146.6159 r?^ + 15.4051 rp - 18.7662 r}%) O^Qiw) 
and 



64217 / 37 6o\ (1) 

[v -v 23 )c ( 8 >{fi w ) 



rf^C^\fi w ) 



9522 



(17) 



39.7055 77*3 -45.4824 7753 + 5.7769 77*3 ) cf\fi w ) 



Explicit expressions for the relevant b — > sj and b — > sg 
matching conditions 



fe=0 



(M, (18) 



4?r 



can be found in Ref. [lj} ■ 

Setting a s (M z ) = 0.118, m t = m t (m t ) = 168.5 GeV, 
fi w = M w = 80.425 GeV and fib = 4.8 GeV, one obtains 



A^Cr(fib 



a s (fib) Y 
4tt / 



and 



17.2 w 0.0051, (19) 



6.4 w 0.0019. (20) 



The correction in Eq. I|19|) causes a suppression of the 
_B — > X s j branching ratio by around 3%. However, 
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one should bear in mind that the size of the correction 
depends strongly on the exact value chosen for jx w . Only 
the total values of C,- (//;,) are guaranteed to become less 
/x w -dependent once more orders in their perturbative ex- 
pansion are being included. At the moment, no complete 
0{a 2 s ) expression for Cj{nb) is available because one 
and only one element on the right hand side of Eq. I|13|) 
remains unknown at this order, namely the element 
^72(m&,Mw) of the evolution matrix that corresponds to 
the current-current operator Q2 = {s L r )^c L ){c L 'y iJ 'b L ). 
It is going to become available once the calculation of 
the four-loop mixing of the relevant four-quark into the 
dipole operators is completed . 

To conclude: We have evaluated the complete three- 
loop 0{cP s ) mixing among the dipole operators, that is, 
in the whole dimension-five sector of the effective theory 
that describes processes occurring much below the 
electroweak scale. Our results are particularly relevant 
for the NNLO analysis of radiative B decays in the SM 
and many of its extensions. 
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